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Abstract 



' We consider an arbitrary translationally invariant chain model with nearest neighbors interaction 

^ ' and satisfying periodic boundary condition. The approach developed here allows a thermodynamic 

description of the chain model directly in terms of grand potential per site. This thermodynamic 
function is derived from an auxiliary function constructed only from open connected sub-chains. In 
order to exemplify its application and how this approach works we consider the Heisenberg XXZ 
model. We obtain the coefficients of the high temperature expansion of the free energy per site of 
the model up to third order. 



o 
o 

1 Introduction 

Chain models, such as the spini Heisenberg models XYZ and XXZ|j], ||] and their extensions j|. U , the one- 
dimensional Hubbard model|| etc., have attracted intense interest due to their property of integrability. 
We also have other models like the one-dimensional d-p model || and its generalizations J?], ||, employed 
to simulate a Cu-0 linear chain with strong Coulomb repulsion. These models are also relevant because 
5_i ■ they may explain high temperature superconductivity. 

The solution of the XXZ model|}| and the one-dimensional Hubbard model|^] at T = are examples of 
the success of the Bethe ansatz. In 1972 Takahashi and Suzuki|n| obtained the exact thermodynamics of 
the XXZ model. More recently the same results were re-obtained by Kuniba et al. jnj using the quantum 
transfer matrix approach. The exact thermodynamics of the one-dimensional Hubbard model was also 
derived by this same method Jl2, [l^, [id). In all those models the thermodynamic functions are solutions of 
integral coupled equations, and it is possible to obtain their behavior in the very high temperature limit. 
However, it is not an easy task to derive a systematic temperature expansion for the thermodynamic 
functions in such temperature range. 

The high temperature limit of any model can be obtained by a standard expansion in terms of 
the inverse of temperature /?|l5|, P = w> wnere & i s the Boltzmann constant and T is the absolute 
temperature. Sometimes, perturbation theory can be applied to obtain the coefficients of the /3-expansion 
method; this approach is also called the high temperature expansion and it has been applied to the t-J 



modelfllfj and to the Hubbard model]!?; 18 ; 191] . In all standard high temperature expansions we always 
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need to calculate the weight of each graph in the grand canonical partition function as well as to determine 
all graphs that contribute at a given /3-order term. 

Recently, we applied the basic properties of the Grassmann generators to calculate analytically the co- 
efficients of the /3-expansion of the grand canonical partition function Z{(3, fx) of self-interacting fermionic 
fields in any dimension ]|(J. We applied it to the one-dimensional Hubbard model^l], We have two 
drawbacks in our method: a large number of graphs must be calculated, as well as to count the number 
of times each graph contributes to each /3-term in the expansion of 2(P,/j,). 

In this paper we restrict ourselves to the study of translationally invariant chain models with first 
neighbors interactions and subject to periodic boundary condition. Instead of calculating the grand 
canonical partition function, we derive a /3-expansion for the grand potential per site, that is an intensive 
quantity. The coefficient of order /3" is exactly obtained for arbitrary values of n. The proposed approach 
is not restricted to self-interacting fermionic models. 

In section 2, we expand the grand canonical partition function as a Taylor series around (3 = 0. In fact, 
the idea behind this method was previously developed, up to order /3 5 , in reference p(| for fermions only. 
Here we generalize the expansion for fermionic and bosonic fields to arbitrary orders of /3" . In performing 
this expansion we write 2(/3, /i) as a sum of traces whose coefficients are known. We obtain the expansion 
of the grand potential per site, and verify that it is an intensive quantity. In this calculation we do not use 
the traditional method of quantum transfer matrix. As an illustrative example of the method developed 
here, in section 3 we apply it to spin-i Heisenberg XXZ model. We calculate the /3-expansion of the 
Helmholtz free energy of the model up to order /3 3 . Taking A = h = we recover the /3-expansion of 
this thermodynamic function of the free fermion model |23]| whereas when t = we recover the limiting 
Ising model p^| . In section 4 we draw our main conclusions. In Appendix A we explain our graphic 
representation of the normalized traces; in Appendix B we extend our results for (K n ) with arbitrary 
n; in Appendix C we write the grand potential per site in terms of an auxiliary function 95(A) that is 
constructed only from open connected sub-chains; and finally in Appendix D we write the function K^"^ 
as sums of normalized traces for n up to 4. 

2 Grand Canonical Partition Function for a Chain Model 

Let us consider a one-dimensional regular lattice (a periodic chain) with N sites, so that the Hilbert space 

N 

of the chain model is simply W- N > = ®TL, H being the irreducible representation at one site, including 
all its degrees to freedom. The dimension of this Hilbert space is dimTi^ = trjv(l). The notation tr/v 
means the trace over all N sites and their internal degrees of freedom, such e.g. spin. 

The grand canonical partition function of a quantum system in the chain with N sites is given by 



where K = H — fiN, with /i being the chemical potential and N being an operator that acts on HS N ^ and 
that commutes with the Hamiltonian of the system. The expansion of 2jv(/3, ji) around /3 = is 



Z N (/3 >f j,)=tr N (e-^), 



(1) 




OO 



(2) 



n=l 
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Let A be any operator that acts on H^ M ^ where M ^ N . We define (A) = j^jftj , for any dimension 
ofH {M) . From now on, we call (A) the normalized trace of operator A. The dimension of the subspace 
7^( M ) is determined by the operator A. 

Using the definition of normalized trace, eq.(||) becomes 



Z N (/3,ti)=tT N (l){l + f^(-Pr { ^}. (3) 

n=l 

Along this section we consider a general Hamiltonian H subject to two constraints: the interaction is 
only between first neighbors and the Hamiltonian H is invariant under translation along the chain. The 
most general operator K satisfying both conditions is 

N 

K = ]TK M+1 , (4) 



where Ki <+i 6 TiS N \ Each operator Kj^+i is defined 



as 



Ki, 2 = Ki, 2 <g> I3 ® • • • <8> ljv 
K 2 , 3 = li ® K 2 , 3 ® I4 ® • • • ® In 

K m+ i = li ig) • • • ® ® K^j+i ® l i+2 ® • • • ® ljv (5) 

and Kj ( j_|_i G 7Yi ®7-^ + i. We use the notation 1^ £ for the identity matrix on the irreducible sub-space 
of the i — th particle. 

Our aim is to get the coefficients (K n ), on the r.h.s. of the eq.(||), in terms of the normalized traces 
of products of operators Kj^+x- Let us calculate explicitly the first three coefficients (n = 1, n = 2, and 
n = 3) of the expansion (||). They will help us to construct the coefficient (K n ) for arbitrary n. We begin 
with n = 1. 

Since the operator Kj^+i acts only on the sites i and i + 1, we have 

tr W (K, !l+ i) = tr 2 (K^j+i ) (tri (1)) N ~ 2 _ ( 6 ) 

We stress out that the traces of the r.h.s. of eq.@ are calculated on the subspace TiS 1 ^ and whereas 
the trace on the l.h.s. is calculated on the complete space ft( N \ Dividing both sides of eq.(||) by trjv(H), 

tr N (K hl+1 ) _ tr 2 (K M+ i) _ 

tr w (l) tr 2 (l®l) " (Km+i) - (7) 

Due to the property of the Hamiltonian being translationally invariant we have: (Kj^+i) = (K^j+i) = 
(Ki j2 ). Taking into account the periodic boundary condition (Kjv.jv+i = Kjv.i), we have that (K) is equal 
to 

N 

(K) = ^(K M+1 ) = 7V(K 1 , 2 ). (8) 

i=l 
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For n = 2 on the r.h.s. of expansion (||) we have to calculate the normalized traces of (K 2 ). From the 
definition of operator K (see eq.(^)), we have 

N 

<K 2 > = {Ki, i+1 K jJ+1 ). (9) 

In the normalized trace (Ki^+iKj^+j.) we have three different cases ( in Appendix A we explain our 
graphic representation): 

i) i =j: 



(K^ +1 ) = (Kf >2 ) — > i i (10a) 

i—1 i i + 1 i + 2 

This equality comes from the invariance under translation and the periodic boundary condition 
satisfied by the model. In the double sum (||) we have N terms of this type. 

ii) i = j ± 1 : 



i =.? + l: (K M+ iK 4+M+2 ) — » i i (10b) 

i i + 1 i + 2 i + 3 



i-i 



(K 



i+l,j+2Ki,i+l 



l) 



i + 1 i + 2 i + 3 



(10c) 



In the double sum on the r.h.s. of eq. (^|) we have N terms of the type of one of these two previous con- 
figurations. Using the cyclic property of traces and translation invariance, we get that the normal- 
ized traces (10b) and (10c) are equal. Therefore, we can write: (K^+iKi-i i) = (K^j+iKj+i^-^). 



Hi) i ^ j and i j ± 1 : 



(K M+1 K,., +1 ) — > ... I i.... (lOd) 

i i + i j j + 1 



The normalized trace (lOd) can be written as (Ki.i + i)(Kj J+ i) = (Ki^) 2 . The number of terms in 
the double sum in eq.(||) that satisfies configuration ( |10d[ ) is equal to N(N — 3). 

The terms that contribute to (IK 2 ), with respective of number of configurations (weight) in the sum on 
the r.h.s. of eq.(^) are 

(K 2 ) = N(K 2 h2 ) + 2^(K l!2 K 2 , 3 ) + jV(^V - 3)(K 1 , 2 K 3 , 4 ). (11) 

The last term on the r.h.s of ( pr| ) can be written as product of normalized traces with n = 1. 

The term (K 3 ) have a richer structure of traces than the two previous ones. Its explicit study will 
help us to generalize the results for arbitrary n. The expansion of normalized traces in (K 3 ) is 

N 

(K 3 )= J2 (K M+ iK 3J+1 K fc , fc+1 ). (12) 
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Similarly to (K) and (K. 2 ), we determine the types of configurations we have in the sum on the r.h.s 
of eq.([l2"l) and their respectively weight. For (K 3 ) we have 

r (K 2 2 K 23 ) (Ki 2 K|,)l 
(K } =N(Ky + hL 2 - 3/ + v hL 2 - 3/ | + 3JV{(K lia K ai 3K3 l4 > + (K 1>2 K 3 ,4K 2 , 3 >} + 

r (K 2 2 K 34 ) (K 12 K 2 4 )-| 
+ W(N - 3){ ^ 12 2l 3 ' 4/ + V h2 2] 3 ' 4/ } + 3N(N - 4){(K 1)2 K 2i3 K4,5> + 

+ (K 1 , 2 K 3 ,4K 4 ,5)} + N(N - 4)(JV - 5)(K 1)2 K M K5,6>. (13) 



From the expressions (g) , ( |ll[) and (|13|) we have that each type of normalized trace (in the expression of 
(K 3 ) we collect them inside of braces) corresponds to all normalized traces that we can write for that set 
of operators. For example, the term: 3(Ki )2 K 2i 3K3 ) 4) + ^{K-i ,2^-3, 4^-2, 3) includes all normalized traces 
that we can write for the set of operators {Ki l2 , K 2) 3 , 1^3,4}. The factors 3 comes from the fact that 
some of them are equals by the cyclic property of the trace. 

To make our notation simpler and generalize our results for arbitrary n, we define the g-trace: 

<K- 1+1 K- a+1 . . . K- m+1 ) g = " l! " 2! n ,'- nm! EWAW • K ^ 2+ i • • • • ' K W TO+1 )>- ("J 

v 

where X)"=i n i = n with m 7^ and the indices ik, k = l..m are distinct among themselves. By definition, 
(V(Ki 1 .i 1 +i , Ki 2i i 2+ i , . . . , Kj mj j m _|_i)) represents all distinct permutations that we can write the n opera- 
tors {Kj^+i , Ki 2 .i 2+ i , . . . , Kj mi j m _|_x}- In the particular case when the operators {K^^+i , Ki 2 .i 2+ i , . . . 
. . . ,Kj mj j m +i} commute among themselves, (K" i 1 ii+1 K" 2 2 i2+1 . . . K™"^ +1 ) g reduces to the normalized 
trace (K- i+1 K- 2+1 ...K^ m+1 ). 

The normalized traces (K), (K 2 ) and (K 3 ) written in terms of the g-traces become 

(K) =iV(K li2 ) g! (15a) 



> =N (K^)g + _/v(K 1)2 K 2i3 ) g + JV(i ! | 3) (K 1 , 2 K 3 , 4 ) g , (15b) 



2! 2! N ' ' * 2! 

and 

(K 3 ) ^ + „ { + (K^U } + ^ (Ki2K23K34)g + 

N(N - 3) r- (K^ 2 K 3 , 4 ) g (K 1:2 K 2 i4 ) g l jV(JV-4) f/iy 
+ 1 2! 2! J + 2! {^i^a.s^.sjg + 

+ (K 1 , 2 K 3 ,4K4, 5 )g} + N ^ N ~ gj (JV ~ 5) (K l|2 K 3 ,4K 5|6 ) g . (15c) 

Some of the g-traces at order n can be written as product of the normalized traces of lower order. 
The latter are already calculated, which means that we have a smaller number of new normalized traces 
to be calculated at order n. To distinguish the normalized traces that appear for the first time at order 
n from the others we use the notation K^"™, where n is the power ((IK™)) of the term in expansion (||) 
to which the g-traces contribute, r is the number of products of the traces that the original trace can 

fn) 

be broken and m is the number of distinct Kj ^+i operators that appear in the g-trace. The terms K^,™ 
with r = 1 are the g-traces that appear in the expansion for the first time at order f3 n . 
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In) 

We begin by defining ~K\ ^ as 



1,1 = n j . ( 16a ) 
K („) = ^ <K?^ 

{rai} 



K l,m = >_> j i j ' ( 16c 

^ ni! n 2 ! ... n m \ 

{m} 

n n 

From now on we use the notation to mean the restriction: ^rii = n and n, ^ for i = 1,2, ..,m. 



The index m satisfies the condition: 1 $J m min(n, iV). In eq.(16c) we are assuming m ^ n. However, 
if m > n we define = 0. 

Since in the definition of the g-traces are included all possible permutation of the n operators, we 
have 

{«;} »=1 * 

For r = 2. the natural definition is 

.. /xrni \ /xrnz \ n .. 



K-(n) _ V^" ^ K 1.2 K 3,4)g _ ( K l,2)g ( K 3,4)g _ V^'V^M" 2 " 



K (") 
^2,3 



ni!n 2 ! f— ni! n 2 ! , , 

{n«} {ni} {n ; } 

/T^ n i T^"2 1<^ n 3 \ IV n ~t T^™3 \ 

\JV! 2 iV 2 3JV4 s^g V- rv l,2 JV 3,4 iv 4,5/e 



(18a) 



, ni! n 2 ! n 3 ! ml n 2 ! n 3 ! 

{«.) (M 



„ ^i,2^ 23 ^4, 5 ; S + ^i^3 4^4, 5 ; g \ = y" ( K M K M +K ^K^) , (18b) 



where we took into account the invariance of the model under translations and the fact that 

when i ^ 7 ± 1, , „ 

\K iti+1 ,K jij+1 \ = < ^ J ' (19) 

I ^ when £ = j ± 1 . 



To obtain the second equality on the r.h.s. of eq.(18b) we rearrange the indices in the double sum. We 



see that these g-traces can always be writing as product of g-traces of lower order which give us 

n m—1 n m 

*S = E" E <Mti - E" E" K SM£ ■ (20) 

{m} 1=1 {m} {rrij} 

The indices rii satisfy the condition n\ + n 2 = n while mi + m 2 = m and 2 ^ m ^ min(n, AT). We are 
assuming m ^ n, otherwise R 2 ™ m = 0. For K33 and R34 we have 

n 

Kg = E'^M^M^i" 3 ' > ( 21a ) 

n 

Kg = ]T" (r^R^R^ + Kf^K^K^ + K^K^K^ ) • (21b) 
{n*} 
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It follows for arbitrary m, where 3 ^ m ^ min(n, N), 



{m} {rrii} 



(n) 

\J 1U1 1 1 L /t. 

We generalize our definition of K^m 



and Kg ^ = for m > n 



K r?m = H 51 II K l"i . (23) 
{n;} {nii} j=l 

where {rii} = {ni,ri2, . . . ,n r } and {mj = {mi,m2, . . . , m r }. Eg. (p3|) is valid for r > 1 and it makes 
possible to write k["™ as products of (see eq.(16c)). The indices r and m satisfy the conditions 

1 r ^ min(n, N) and r ^ m ^ min(n, N). We are assuming n ^ r and m, otherwise kJ."!, = 0. 
In terms of K^ m , the normalized traces (IK), (IK 2 ) and (K 3 ) are simply 

(K) = NK^l , (24a) 
(5!> = ™g + ™g + ffi^K», (24b) 



<|> . W K™ + W K™ + WK j + S^K g + ffipiiKg + ffi^iK^iKg . ,24c, 

In the Appendix B we write (K 4 ) in terms of K^i and the coefficients of K^and Kg™^ for 

n = 1, 2 and 3 are rewritten in terms of binomial coefficients. In this appendix we also extend the previous 
results for arbitrary n. 



For arbitrary n, eq.( |B.8| ) gives the expression of (K n ), that is 



r—1 m—r 



The notation [n, N] means the min(rt, N). Differently from our previous work|2(| [2l], [22) , in cq. (p5|) wc 
already have the weight of each set of sub-chain Kf-^j in (K n ). The coefficient that multiplies K^L in 
eq.(|2^) is independent of the particular Hamiltonian under consideration. 

For n < N the trace will be calculated on a maximum subspace given by HS n+1 \ whereas for n ^ N 
due to the periodic boundary condition the operators act on whole Hilbert space TL^ N \ therefore we can 
say that the operators act more than one period in the periodic chain. 

Rewriting the coefficient of K^i in eq. ( p5| ) in a more convenient way, 

r=l m=r k=l ^ / \ / 



EC 

k=l v 



&k,n (26) 



with 



[n,N] [n,N] . 

vnrq m+r i mkw. (27) 



r—k m—r 
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From eq. ( |27|) we have that the biggest possible value of k is mm(n, N) while there is no restriction on 
n. The function &k,n has the property, 

n fe 
{».) *=1 



where {n^} = {rti, ri2, ■ . ■ , n^}. Replacing eq.(|2q) in eq. (j3[) we get 

JV / , ,s oo 

^05 lA i)=tr J ,(i ) {i + x;( JE(-^"M- (29) 



fe— 1 n— fe 



From property fl28[ ) and recombining the summations conveniently, eq. (|29|) is finally written for any value 
of N as 



Z N ^,fi) = {tv 1 (l)(l + J2(-(3) n ^,n)} ■ (30) 

n=l 

The grand potential per site W(/3, //) is 

W JV (AM) = -^ln^(AM) = -^ln(tr 1 (l)(l + ^(-/3)^ 1 ,„)). (31) 

n— 1 

This expansion can be used to obtain analytical results of chain models in the high temperature limit. 

We can be misleading by the notation R\ t n in eq.(|3l|) since in it are included not only connected 
sub-chains but also disconnected ones. For finite N, eq.(|3l]) is the simplest form of the grand potential 
per site since in this case n can be bigger than TV which means to go around the chain more than once. 

2.1 The Thermodynamic Limit 

Our main interest is to get the thermodynamical behavior of physical quantities that characterize physical 
systems. Those are attained only in the thermodynamic limit (N — > oo). We first take the thermodynamic 
limit and after makes the expansion (^|), so that we never goes along the chain more than once when we 
calculate the normalized traces. In this thermodynamic limit the r.h.s. of eq.(^0|) is equal to the highest 
eigenvalue obtained from the associated quantum transfer matrixp4|] to the quantum system of interest. 
From eq.(|T]), the grand potential per site W((3,n) in this limit is 

WGM = - \hn—\nZ N ([3^) = — In (tn(l)(l + £ (-/?)"*!,„)) . (32) 

n— 1 

Eqs.(|3T|) and (]32j), look alike. Their difference relies on the fact that in cq.(]3^) the biggest sub-chain that 
contribute to &\^ n has at most the length of the chain. Eq.(|32|) confirms that the grand potential is an 
intensive quantity. 

Certainly our aim is to be able to calculate the smallest possible number of open sub-chains to 
obtain the grand potential per site. But the definition of (see eq.(27)) includes connected as well 
disconnected sub-chains. In eq. (|32|) we have the following sum to calculate 

oo 

£ = £(-/3)»S lfn . (33) 

n=l 
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Let us define 



oo 

£ 



(-/3)"K^ 



(34) 



as a summation of powers series of 0. The expression of is given by eq.(16c). We point out that 

only connected open sub-chains contribute to the functions T m . We define the function <p(\) in terms of 
the functions T rn as follows 



<pW = E 



A' 



(35) 



where A is a parameter. In [Appendix Cj we show that the 

>(A)" +1 



d" 



oo 

E dA 

n=0 



(n+1)! 



(36) 



A=l 



From eqs.(p7|) and ( |32[) we have that the connected and disconnected open sub-chains contribute to 
the grand potential per site. In the thermodynamic limit we have from eqs.(^4|)-(^6|) that all these sub- 



chains can be derived only from a sum of connected open sub-chains (for details see Appendix C ). We 
finally write the grand potential per site simply as 



W(/3,M) = ~{ln(tn(l))+ln(l + 0}. 



(37) 



Due to the fact that £ can be derived from function ip(X) one has the consequence that a much smaller 
number of terms has to be calculated at each /3-order in /3-expansion of W((3,fi). The weight of each 
g-trace in the eq.(|37|) is already included in the definition of £. 



3 The Heisenberg XXZ Model 

Let us consider the Hamiltonian of the well known anisotropic one-dimensional Heisenberg XXZ model 
with spin- ig§ 

1 N 

3=1 

where h is the external magnetic field and (crj, crj, cr|) are the Pauli matrices at the jth site on a periodic 
chain with N space sites and A is called the anisotropy parameter. For A > (A < 0) we have a repulsive 
(attractive) interaction core. The case A = 1 (A = — 1) corresponds to the isotropic antiferromangetic 
(ferromagnetic) Heisenberg model with fully polarized ground state. 

Through the Jordan- Wigner transformation |^5j , the Hamiltonian ( |38|) is mapped on a spinless fermionic 
model, whose Hamiltonian is 

N A 
H = (*(o]-°i+i + a ]+i a j) + Vn J n i+i + En i) + N ( h + -J )' ( 39 ) 

3 = 1 
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with V = 2 A, E = —2h — 2A, where fii = ojaj and a] (Oj) is the fermionic creation(annihilation) operator 
at site j. The operators a\ and Oi satisfy the usual anti-commutation relations of the fermionic fields. 
The hopping constant t is included in Hamiltonian ([59]) only to help us counting the powers of the terms 
that contribute to a given order j3 n in the /3-expansion of the Helmholtz free energy per site. Throughout 
our calculations, we take t = 1. The term N(h + —) leads to a shift in the energy of the ground state of 
the model. 



3.1 High Temperature Expansion 



In this section we calculate the /3-expansion of the Helmholtz free energy per site of the XXZ model, 
whose Hamiltonian is given by eq. (|39|) , up to order /3 3 . 

In order to apply the method derived in section 2, we write Hamiltonian (f38|) as: H = J2 i=1 H i]i+ i, 
where 



H M+1 = E M+1 + T+ +1 + T7. +1 + V M+1 , (40) 

^i.i+l = "*»"*+!) A i,i+1 



and each term on the r.h.s. of cq.(|40|) is defined as: E^j+i = \E{rii +nj + i), T- i+1 = ta]a.; + i, T 



ta\ +1 ai and V M+1 = Vnin l+1 . 

Our aim is to obtain the Helmholtz free energy per site of the XXZ model in the thermodynamic 
limit. For doing that we need the auxiliary function <p(X) (see eq.(p9)) of the XXZ model. This function 
is calculated from the functions T m that are written in terms of (see eq.(Q)). Since we intend to 

calculate the /3-expansion of the Helmholtz free energy per site of the XXZ model up to order /3 3 , we 
need to obtain the analytical expressions of T±, . . . ,T±, namely, 

I\ = -I3K.Q + /3 2 K^ - /3 3 K^ + /3 4 K^ + G(/3 5 ), (41a) 

T 2 = /3 2 K^ - /3 3 K^ + /3 4 K^ + 0((3% (41b) 

r 3 = -/3 3 K^ + /3 4 K$ + 0(/3 5 ), (41c) 

r 4 = /? 4 K( 4 l + 0(/3 5 ). (41d) 

In appendix D the functions i n eqs.(^ll) are written as sums of normalized traces. 

In particular, for Hamiltonian (|39| ) the hopping terms (T^" i+1 and T~ i+1 ) in the normalized traces 
only contribute to those functions K^, where we have H' i+1 with I > 1. Since we are calculating traces, 
the number of T^ +1 's in a given normalized trace has to be equal to the number of terms T~ i+1 's in it. 
Otherwise, the normalized trace is null. 

Calculating explicitly the normalized traces that contribute to ( see appendix D) for the XXZ 

model, we obtain 



• m = 1 



K^ = y + T , (42a) 

tAA) i 4 EH 2 3£ 4 VE 3 E 2 V 2 EV 3 V i , An „ 

K ] = 1 1 1 h 1 h — . 42d 

1,1 48 32 256 24 16 24 96 ' 



E V 
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(2) 
1,2 

(3) 
1,2 

(4) 
1,2 



• to 



K 
K 

• to = 4 



(3) 
1,3 

(4) 
1.3 



5E 2 SEV 
"T6~ + 8 



V 



4 

24 " 



3E 3 
IF 



V 3 Et 2 17VE 2 7EV 2 Vt 2 
5E 2 t 2 95£ 4 



24 



16 



+ 



32 16 8 

25V E 3 33E 2 V 2 



768 



64 



+ 



64 



5EV 3 
16 



7^4 ^2y2 

96 + 24 



23VE 2 EV 2 



64 



41£ 2 t 2 59E 4 85VE 3 101E 2 V 2 7EV 3 7t 2 V 2 23EVt 2 



192 



256 



128 



128 



16 



96 



96 



K 



(4) _ 
1,4 — 



29.E 4 19VE 3 5E 2 V 2 5EV 3 



256 



64 



16 



32 



Yt 

32 ■ 



EVt 2 



3V 4 
"32"' 



(43a) 
(43b) 
(43c) 

(44a) 
(44b) 

(45a) 



Substituting those expressions in eqs.(f4l|) (taking t = 1) and those in eq.(p5|) we obtain the auxiliary 
function <p(X) for the XXZ model. Using the algebraic language MAPLE 5.1 we calculate the Hclmholtz 
free energy per site of the model up to order (3 3 , that is, 



W) = ^ + Wi(0) + W F (J3) - j/3 2 + \(h 2 + ^)/3 3 + 0(/? 4 ), (46) 
where Wj(/3) is the Helmholtz free energy per site of the Ising model[^4|, |2(| up to order /3 3 , that is, 

and Wf(/3) is the Helmholtz free energy per site of the free fermion modcl[^6|, |23| up to same /3-order, 

W,G9) = -tp- *-0> + ((3 4 ). (48) 

We rewrite the Helmholtz free energy per site (eq.([l6|)) using the constants of Hamiltonian in reference 
p7f . Comparing our results with eq.(6.5) of this reference, we conclude: i) both results agree up to order 
j3; ii) our coefficient at order (3 2 is: — J 3 cos(7) + Jcos(7)/i 2 which means that we have two misprints in 
eq.(6.5) of reference |p7| , that is: the sign of — J 3 cos(j)f3 2 and is missing a power of J multiplying h 2 
(the two terms that contribute to order /3 2 in eq.(6.5) have different dimensions); in) it is missing the 
term — J 2 cos 2 (j)h 2 f3 3 in eq.(6.5) of reference |27j (we point out that without this term we do not recover 
the limiting case of Ising model |m| ^6| (t = 0) from XXZ model), morever the other terms at order /3 3 
in eq.(6.5) have the following misprints: — j,J 4 j3 3 should be \ J A 3 and the term — | J 4 cos 2 (7)/3 3 should 
be | J 4 cos 2 (7)/3 3 . With the previous corrections we get the Helmholtz free energy per site of the Ising 
model and the free fermion model as limiting cases of the XXZ model. 

In reference [27 Destri and de Vega obtained the Hclmholtz free energy per site of the XXZ model as 
a non-linear integral equation (NLIE), as similarly obtained by Klumperp3|. Once our disagreement with 
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the result of reference |27j is proportional h 2 at order /3 3 , we calculate the /3-expansion of the magnetic 
susceptibility x(0) (x(P) — — 9 obtained from eq.(|4rjf). In figure 1 we plot the curve of the magnetic 

susceptibility derived from NLIE in reference |23|, ^7), the /3-expansion of x(P) derived in reference f27f 
and our own /3-expansion of x{P)- From figure 1 we see that our result with A = \ and A = ^ coincides 
with the one derived from the NLIE in reference |2j| |27j] for up to 0.3, while de /3-expansion of Destri 
and de Vega does not. 



Magnetic Susceptibility Magnetic Susceptibility 




(3 



Figure 1: The pointed line corresponds the /3-expansion ([l6|), the dashed line represents the /3-expansion 
obtained by Destri e de Vega[^7) and the solid line corresponds the numerical solution of the NLIE in 
reference [23, 271]. 



As a final check of our /3-expansion (|46|) we plot in figure 2 the third partial derivative with respect 
to (3 of the function /3W(/3). We compare the numerical solution for this function obtained from the 
NLIEg |2p, in the same reference, and our results derived from eq. (p6|) . In the region /3 « 0, the 
function ^3-(/3W(/3)) is a straight line. From figure 2 we see that the result derived from eq.(^6|) is 
identical to the one from the NLIE for /3 up to 0.1. For h ^ in figure 2 we get that the function 
^p- (/3W(/3)) derived from the /3-expansion in reference does not even touch the exact curve. 
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p p 

Figure 2: We plot the function ^(/3W(/3)) versus p. The pointed line (a straight line) corresponds to 
the results obtained from (|46|), the dashed line the /3-expansion obtained by Destri e de Vegaj27j and the 
solid line the numerical solution of the NLIE|2^, 

4 Conclusions 

There are several interesting one-dimensional models with translational invariance, with interactions 
between first neighbors only, satisfying periodic boundary conditions. They are generically called chain 
models, some which having the property of integrability and to which the Bethe ansatz has been applied 
to solve them exactly at T — and at finite T. Their thermodynamical properties are obtained through 
coupled integral equations which demand numerical analysis to be solved (as examples see references 

In this paper we present a new analytical method to get the /3-expansion coefficients of the grand 
potential per site, for translationally invariant chain models with nearest-neighbors interactions. Differ- 
ently from our previous works [f^, [2l], ||^| we calculate directly the /3-expansion of the grand potential 
per site W(/3, fi). The weight of each sub-chain in W(/3, /x) is obtained, and we show explicitly that the 
grand potential per site is an intensive quantity. In the thermodynamic limit (N — > oo) we show that 
W(P, /x) can be derived from an auxiliary function ip(X) that is written only in terms of open connected 
sub-chains. The present approach gives a /3-expansion of W(/3, fi), whose coefficient of order P n can be 
obtained exactly for arbitrary value of n. The existence of this auxiliary function allow us to get higher 
order terms in the /3-expansion of W(/3, fi) than we were able before, in references |^0|, [2l], ^]. The coeffi- 
cient of order /3 n is analytically obtained, and it is exact. The number of terms (traces) to be calculated 
when derivig each coefficient in /3-expansion of W(/3, fi) is much smaller than before. 
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The present method applies to chains in both real space and momentum space, and the result given 
by eq.(|37|) is valid for both classic and quantum models (bosonic and/or fermionic fields). 

As an illustrative application of the method, we consider the spin-i Heisenberg XXZ chain model and 
derive its Helmholtz free energy per site up to order /3 3 . Our /^-expansion (^) recovers the Helmholtz free 
energy per site in two limiting cases: the Ising model (t = 0) and the free fermion model (A = h = 0). 
We compare our high temperature expansion of the W(/3) eq.([46|) of the XXZ model with the one derived 
by Destri and de Vegap7|]. We show that this last one has a misprint in eq.(6.5). 

As a final word we want to point out that the present approach can be well applied to models that 
do not have the property of integrability. Besides, once we have closed expressions for the coefficients of 
the /3-expansion of W(/3) the method can be implemented through a procedure in an algebraic language 
to get high /3-order terms of the high temperature expansion of W(/3) for the model under interest. We 
hope to report in the future progress in this direction. 
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Appendix A Graphic representation of normalized traces 

We have a better visualization of the normalized traces that contribute to (K n ) through a graphic rep- 
resentation. In the graph, the operator Kj^+i is represented by a straight line that links the sites i and 
i + In the case we have a product of operators K^f+i's for each operator a line is drawn at a different 
step of the ladder. It is important in the normalized trace the order in which the operators K^j+i appear 
in the string. We take this order into account in our graph representation of the normalized traces. To 
exemplify how to take into account the order of the operators in the string, we consider the following 
examples: 

(K M+1 ) — » . . (A.l) 

i - 1 i i + 1 i + 2 



i i + 1 i + 2 i + 3 



(A.2) 



(K^i+iKj_ li K i+li+ 2) 



1 i i+1 i+2 



(A.3) 



(K,:,t+iK i+1/j+ 2 



K; 



i) 



i+1 i+2 



(A.4) 



The graphs (A.1)-(A.4) are connected graphs, which means that the normalized trace can not be written 
as a product of normalized traces of smaller number of operators K» i+i. 
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Appendix B Calculation of (W 1 ) 

In this appendix we give the explicit definition of K^ 4 m in terms of which we write (K 4 ) . The expression 
of (K 4 ) plus eqs.(E4f) will allow us to write down (IK™) for arbitrary n. 



We define k[ ^ m as 



K$3 = ^, (B.l„i 



T ,(4) ( K l,2 K 2,3)g ( K l,2K 23 ) g (K 12 K 23 ) g , D1U 

K - = 3! + 3! + 2! 2! ' (B ' lb) 

(4) (Kj |2 K 2 , 3 K 3 , 4 ) B (K 1 , 2 Kl i3 K 3 ,4) g ^ (K 1;2 K 2 , 3 K^ 4 ) g 

KU = ^ + + g , (B.lc) 

Kg - (K 1 , 2 K 2i3 K 3 , 4 K 4:5 )g. (B.ld) 



The other traces can be written as products of traces of lower order and they are 
(K x 2 K 2)3 K 4 ,s)g (Ki !2 K 2 3 K 4i5 )g (Ki )2 K 2 , 3 K 4 5 ), 



T -ii! \ K i,2 K 3,4) g (Ki i2 K 34 ) g (K 12 K 34 ) g 

K, , = ; — - - — ■ — - J — , (B.2a) 

2,2 3! 3! 2! 2! v ' 



K (4) 
^2,3 



2! 2! 2! 

(K? 3 K 3 ,4K4, 5 }g (K ll2 K§ 4 K 4 , 5 ) g (K 1)2 K 3) 4K2 5 ) 



2! 2! 2! ' (R2b) 

K$ =(K li2 K 2 , 3 K 3 , 4 K 5 , 6 ) g + (K li2 K 3 , 4 K 4 ,5K 5 , 6 }g + K 1)2 K 2)3 K 4>5 K 5)6 ) g , (B.2c) 
(4) _ (Kf ;2 K 3 , 4 K 5 , 6 ) g (K 1 , 2 K|^K 5 , 6 ) g (K^I^jKjU 

3,3 2! 2! 2! ' ' 

K$ =(K li2 K 2 , 3 K 4 , 5 K 6 , 7 ) g + (K 1 , 2 K 2i3 K 4 , 5 K 6)7 ) g + (K 1)2 K 2)3 K 4i5 K 6)7 ) g , (B.2c) 

Kg =(K 1 , 2 K 3;4 K 5 , 6 K 7 , 8 ) g . (B.2f) 

From those definitions, the normalized trace (K 4 ) is written as 

(^l , ^(4) , ^(4) , ^(4) , N(N-S) jAi) , iV(iV-4) T , (4) , A^V-5) (4) 



4| = JVKft + NK% + NK% + JVKft + " 2 , ' K% + - ^ ' K% + " - { ' K 2A 

N(N-i)(N-5) jAi) N(N-5)(N-6) jAi) N(N - 5)(JV - 6) (TV - 7) (4) 
3! 3 ' 3 3! 3 < 4 4! 4 < 4 ■ 



We use the following notations to write the coefficients of eqs.fcj) and eq.(B.3): 



riJV 






W,2 




^2,2 — 


^1,3 




^2,3 — 


riJV 
u i, 4 


= iV, 


^2,4 — 



N(N- 


3) 


2! 




7V(7V- 


4) 


2! 




7V(7V- 


5) 



2! 



C 



3,4 



7V(7V-4)(7V-5) 
3! 

7V(7V-5)(7V-6) 
3! ! 



C 



N 
4,4 



iV(iV-5)(iV-6)(iV-7) 
4! ' 



(B.4) 
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Replacing the definitions given by eqs^B^) in cqs.(gj) and (BJ3), we obtain a compact way of writing 
the normalized traces 



(K) = C^K« 

TK-2\ 



2! 
3! 
4! 



r N t^ 2 ) _l x^( 2 ) _i_ pJV T^( 2 ) 

^1,1^1,1 T ^1,2-^1,2 ^2,2^2,2 ' 

r N k (3) 4- r N k (3) -i- r N k (3) + r w k (3) + r w K" (3) 

^1,1^1,1 ' ^1,2^1,2 ' ^1,3^-1,3 ' ^2,2^-2,2 ' ^2,3^2,3 



C N K (3) 
^3,3^-3,3 ' 



pJV K (4) , r N K (4) pAf K (4) p AT K (4) r N „ (4) pAT _(4) pAT _(4) 
^1,1^1,1 T ^1,2^1,2 "•" ^1,3^1,3 b l,4 A l,4 + ^2,2^2,2 "+" ^2,3^2,3 + ^2,4^2,4 

, r N K (4) , r N K (4) , r N K (4) 
~r ^3,3^-3,3 "T" ^3,4^3,4 "+" ^4,4^4,4 ■ 



From eqs.(B.4) we see that coefficient that multiply is independent of n and it is equal to 

CZm = N, Km, 

/VY/V — m — 1) 

— , 2 ^ m, 



pA _ 
^2,m — 



2! 



jv _ AT(AT-to - 1)(AT-to-2) 

^3,m 



3 $5 m, 



C 



w _ N(N -m-l)(N -m-2)(N -m-3) 
~ 4! 



4^ m. 



Afterwards for any r ^ to ^ min(n, AT), we obtain a general relation from eqs.(B.6) given by 



N NfN-m-1 
r ' m r r - 1 



fe=i 



such that 



r—1 m—r 

We are using the notation: [n, A 7 ] = min(n, N). 



r - 1 



(B.5a) 
(B.5b) 

(B.5c) 
(B.5d) 

(B.6a) 
(B.6b) 

(B.6c) 

(B.6d) 

(B.7) 



(B. 



Appendix C Relation between the auxiliary function (p(X) and 

In the eq.(|34|) we define the function T m , 



r m = £ (-/3)"k< ; 



(„) 

m ' 



(C.l) 



and the auxiliary function ip(X), 



(C.2) 



m— 1 
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In this appendix we show that the grand potential per site, in the thermodynamic limit, can be obtained 
from derivatives of powers of the function <£>(A). 

To calculate the grand potential per site (see eq.(|32|)) we need to perform the sum 



(C.3) 



Substituting the expression of &i m/n (eq.j27|) with k = 1) in £ and taking into account that in K^ n m we 
have to have n > m, eq.(C.3) becomes 



oo n 



oo n 



7i— 1 m=l n—2 rn — 2 * 

oo 

r=l 

where y r is such that 



oo n 



K 



2,m 



5EEN) 



n— 3 m— 3 



3,m 



, , y r+1 oo n 

^EEH) 



n— r m—r 



m + r — 2 
r - 1 



K (n) 



(C.4) 
(C.5) 

(C.6) 



We drop the upper limit [n, N] condition in the sum over m (see eq.(p7|)) because now we are taking first 
the thermodynamic limit (N — > oo), which means that we always have n < N. 



Applying the equality 



OO 11 



oo oo 



EE — EE 



(C.7) 



n—r m—r 



m—r n—r 



on eq.(|C.6|), we obtain 



(-l^+i « fm + r-2 



e ( r _; ; m-® n K% 



(C.8) 



From eqs.(C.l) and (C.2) we have 



^) = EiE(-^ 



(n) 



(C.9) 



m— 1 n— 1 



Our next step is to take powers of function <p(\). We begin with 99(A) 2 



mi, 777,2 — 1 77-1 ,772 — 1 



n m 



E ^ Eh^E" E" K S K S 

m— 2 n—2 {r, 

00 1 00 

EiE(-^ 



{n;} {mi} 



-in) 



(CIO) 



TTi—2 n—2 
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For arbitrary r we obtain 



and we note that 



OO oo 



^ X m ^ 

m—r 7i— r 



p( A )|a=i = Pi 
d MA) 2 



AW 2! 



A=l 



Y2 



(C.ll) 



(C.12a) 
(C.12b) 



Finally we have 



clA 



- 1 MA) r \ 



A=l 



r=l 

= V- 



</>(A) r 



A=l 



(C.13) 



which shows that the sum ( |C.3| ) can be obtained from a function that contains only open connected 
sub-chains. 



(C.12c) 



Appendix D Useful relations between and the normalized 
traces 



In section 3 we calculate the /3-expansion of the Helmholtz free energy per site of the XXZ model up to 
order /3 3 . In order to do so, we need to write the functions kJ™2i ( see eq.(16c)) in terms of the normalized 
traces instead of the g-traces. 

For any quantum system driven by a Hamiltonian with interactions between first neighbors only, 
invariant under translations along the chain, and null chemical potential, we have: 



• m = 1 



(i) 
1,1 

(2) 
1,1 

(3) 

^ " 3! 

(4) _ 1 
1,1 



(Hl,2 

1 

2! 
1 



<H£ 2 >, 



(H?, 2 ), 
4! < H U- 



(D.la) 
(D.lb) 

(D.lc) 

(D.ld) 
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Kr2 = (Hi,2H 2i 3), 



1,2 2! 

'(4) 



(Hi i2 H 2:3 }+ (Hi^H^g) 



K$ - (H? i2 H 2 ,3)+ -(2(H? i2 Hl j3 }+ <H li2 H 2i3 H li2 H 2 ,3» + (H 1;2 H| 3 ) 



2! 



• m = 3 



= (Hi j2 H2,3H 3j 4), 



(D.2a) 
(D.2b) 

(D.2c) 
(D.3a) 



(4) 



2 ,1 1 

(H 12 H2 j 3H 3i 4) + -(H li2 H2 ;3 H 3i 4) + (H li2 H2 i 3H 3j4 ) + -(H^H^Hs^B^s) . (D.3b) 



• m 



K^l — (Hi j2 H 2i 3H3 j 4H4 )5 ). 
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